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Current-driven domain wall motion in thin ferromagnetic wires 
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The coupling between a current and a Bloch wall is examined in the half-metal limit of the double 
exchange model. The conduction electrons transfer angular momentum to the Bloch wall with 
100% efficiency in the absence of pinning. The wall is displaced without distortion with velocity 
proportional to the current. In the presence of a pinning potential either the angular momentum 
is destroyed by the perpendicular component of the anisotropy field or is converted to coherent 
magnons. A moving wall has its velocity reduced by pinning. The expression for the velocity agrees 
surprisingly well with experiment. 

PACS numbers: 75.60.-d, 73.40.Cg, 73.50.Bk, 75.70.i 



Spintronic devices have great technological promise 
but represent a challenging problem at both an applied 
and fundamental level. It has been shown theoretically [1, 
2] that the direction of a magnetic domain might be 
switched using currents alone. Devices designed to use 
this principle often consist of multilayers of magnetic and 
non-magnetic materials. The advantages of similar de- 
vices based upon the current induced displacement of 
Bloch wall are simplicity and the fact that the switch- 
ing current is smaller [3-5]. Experimentally the current 
induced displacement of a Bloch wall has been clearly 
demonstrated and in a recent experiment the velocity of 
the displaced wall was measured [5]. 

In this Letter is developed a first principles theory for 
the coupling between a current and a Bloch wall based 
upon a standard model for ferromagnets, i.e., within the 
double exchange model. Great care is taken in order 
to ensure that angular momentum is conserved. Using 
the s — d-exchange model, Berger[2] shows the polariza- 
tion of the s-electrons rotates as it passes through a wall. 
The reaction force might then be assumed to cause wall 
movement. However, as will be shown, the situation is 
complicated. Angular momentum given to the wall can 
be taken up either by (i) the displacement of the wall, (ii) 
the production of coherent spin waves or (iii) a rotation 
of the wall so that the angular momentum is absorbed 
by the anisotropy energy. These issues are most eas- 
ily addressed in the half-metal limit. This is technologi- 
cally pertinent since it corresponds to the most efficient 
devices. Most work has been performed on Permalloy 
(Nigilfeig). For this material, the current j has a polar- 
ization p ps 0.7[5] and so the present description based 
upon p = 1 is not hopelessly unrealistic. When mak- 
ing comparison with experiment, the difference will be 
accounted for by adding factors of p as appropriate. 

The standard double exchange model is 



n = -* E ( c W + H - c ) + U J2 n w 

<ij>a i 

J2( J HSi • ^ + AS lz 2 - K ± S iy 2 ) -J J2&- ^ W 



where Si is the localized spin operator and c ia creates 
an electron with spin a at site i. The wire is assumed 
to extend along the z-direction. That this an easy axis 
implies A > 0. 

Corresponding to highly correlated electrons, the limit 
U ^> t is taken. This permits a spinless majority electron 
operator c\ = ctj. to be defined, and in terms of this the 

minority cji = s~c\ where s~ is the conduction electron 
spin operator. A large Hund's rule coupling | Jjj\ 3> t is 
also assumed. With this only the largest total angular 
momentum manifold is relevant. At a site occupied by a 
conduction electron only states with J = S+ 1/2 need to 
be accounted for and it follows that Si w Si/(2S). This 
is not an assumption about local equilibrium of the con- 
duction electron but rather is an immediate and rigorous 
consequence of the Wigner-Eckart theorem. 

The Holstein-Primakoff transformation [6] Si e = S — 
b\bt, S+ = (25 - &l&i) 1/2 &i » (25) 1/2 &, is used to quan- 
tize the spin degrees of freedom. In principle, the axes of 
quantization are arbitrary, however, as the standard the- 
ory of ferromagnetic magnons illustrates, in fact the lo- 
cal axis of quantization must be taken along the classical 
equilibrium direction, i.e., the spin direction in the limit 
S — > oo. When this is the case the approach generates 
and expansion in 1/5. This local equilibrium direction is 
specified by the Euler angles 0* and . In this Letter the 
solution will most often be time dependent. Rather than 
setting up the complicated equations of motion, here ex- 
tensive use will be made of rotating frames. The aim is 
to reduce the problem to that of static equilibrium but 
within the time dependent frame. The time dependent 
solution is then obtained upon reverting to the labora- 
tory frame. 

In the absence of a wall, for large U, the kinetic 

energy term is T — -tJ2(ij)( c l c j + s ~i c \ c j s ^) + H - c - 
The Bloch wall is generated through via the static ro- 
tations by 8, and 0. For this large U regime, the usual 
SU(2) rotation of the conduction electron operators sim- 



plifics to c] -> (cos [0/2] 
tated f w -*E<y)Cos[(0i 



sin[0/2] s, )cj and the ro- 
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motion arises principally via TL t 



4c jS f) 



H.c. Charge 

t£ w> t«ctc J . + ff.c. 
where ijj = fcos[(#j — #j)/2]. The reduction of tij in the 
Bloch wall implies a barrier, however this has a height 
~ SA which is taken to be negligible compared to Ef- 
The solutions of Ht are, to a good approximation, plain 
fc-states independent of the wall position or its motion. 

Using the Wigncr-Eckart theorem, the part of T linear 
in a? reduces to: 



H t j — - 



86; 



2(2Sy/ 2 dz 



aY J 4ci + i{b i -b\ +1 )+H.c, (2) 



which is the key spin-charge current interaction. The, bi- 
linear in sf , part of T is also of importance. When act- 
ing on eigenstates of Ht these lead to a renormalization, 
J = J + (x't/2S 2 ), of the total exchange coupling J. 
The effective concentration x' = (1/N) X^(cos k z a) ng, 
assuming isotopic hopping. 

A small wall translation by Az amounts to a 
small rotation generated by R(Az) = Y\ i cxp(i(9 i — 
8 l+1 ){Az/a)S ly ) w {\ + {l/2)Y. l {de i /dz)Az{2S) 1 ' 2 {b\- 
bA) = 1 + S 1 / 2 (2A/J) 1 / 4 (Az/a){bl - bj which defines 
bl = (J/8A) 1 / 4 J2 l (d9 t /dz)ab\. That the translation 
R(Az) leaves the energy unchanged implies that b^ cre- 
ates a zero energy excitation localized in the wall. This 
special frj^-mode is key to a understanding of not only 
the dynamics but also pinned static solutions of a Bloch 
wall. It is of the nature of a discrete Goldstone boson 
associated with the translational invariance of the wall. 

The effective, Ti w E{O u <fc\ + H t + H s + H tJ , 
where E{9i,<j>i} is the classical energy functional, Tit is 
given above and Ti s is the magnon Hamiltonian. The 
Bloch wall structure can be determined by minimizing 
E{ei,(j>i}. The standard result is 6(z) = 2 cot" 1 e - {z / w ^ 

with = for K± > 0. The wall width w = a (^) 1/2 - 
The local axes change from 9 = to 9 = it with increas- 
ing z and the wall lies in the x — z-plane. 

The conservation of the spin angular momentum J z 
is a central issue. If the coupling to the contacts is ig- 
nored and K± = 0, it follows that [J z ,Ti\ = 0. A wall 
with speed v implies J z — hSv/a and, in fact this must 
be equal to the surface integral J s ■ dA where s is the 
spin current. The conduction electrons have their spin 
reversed in the wall implying a contribution J s- dA oc /, 
the charge current. Such a finite surface term gives rise 
to a local term V • s oc j, the charge current density, in 
the equation for dSi j dt and will lead to wall motion even 
when K± = 0. 

Whatever the value of K±, the translational invari- 
ance of the Hamiltonian implies that the wall will move 
with a constant velocity in the presence of a charge cur- 
rent. In order to prove this, use is made of a rotating 
frame. Consider again the current carrying states which 
are eigenstates of — etc-. The spin-charge coupling 



reduces to 
Titj 



.hja? 
'' 2eS 



£ga(2S)V^ 



4). (3) 



Here {2S) 1 / 2 {l/2i){b l - bj) w Sf y defined to be strictly 
perpendicular to the instantaneous axis of quantization. 
The only derivative which is finite during this short pe- 
riod is the similar instantaneous Sf x — S^g- J2i ^i a - 
This derivative can also be made to be zero by a suitable 
rotating frame, i.e., via r = exp(i(d9i/dt)Sf y t) so that 



rS, 



and Sf, 



S ,, 



i(d0i/dt)[sLsr\ 



Six — i(d9i/dt)S. The wall is therefore stationary, when 
the local axes of quantization rotate according to 



dt 



v - 



d6i 



2eS' 



(4) 



The solution of this is 9a 



2 cot" 



-(z/w) 



8i(z — vot) where 6i(z) 



is the static solution. In the laboratory 
frame, the wall is translated without distortion. The ve- 
locity wo [5, 7] is such that the net conduction electron 
spin current which enters the wall is exactly compensated 
by the similar current implied by its uniform displace- 
ment. 

The simplest model for the pinning of the Bloch wall 
is to reduce by a small amount the anisotropy energy at 
the origin. The added term Ti p = +SASq z 2 generates 
a pinning gap E p . The problem of such a pinned wall 
is surprisingly complicated. Our results developed below 
show, for finite K±, the wall tilts by an angle <f> and the 
critical charge current j c is determined by the condition 
that E p = 0. However, when K± is negligible, the wall 
generates magnons which take away the conduction elec- 
tron angular momentum given to the wall. In this case, 
j c is found to reflect the equality of the rate of creation 
of linear magnon momentum to the maximum force gen- 
erated by the pinning. 

The largest pinning gap E p arises when the wall is cen- 
tered at the origin, i.e., with (9o = 7r/2. Explicitly then, 
Soz — > Sq x and the diagonal part of Tip = 2S5Ab' b . The 
value of Ep is determined by evaluating the commutator 
[[Hp, b w ], 6^], which is equivalent to taking the expecta- 
tion value with respect to the 6|„-magnon wave function, 
and, with j = 0, gives, 



E p — S 



2A 

T 



1/2 



SA. 



(5) 



The wall can now tilt so that <f> is finite. The anisotropy 
K±J2i^iy 2 contains — 3SK± sin 2 9i sin 2 4>ib\b i which 
gives the fe^-mode an energy shift of —2SK± sin 2 <f>. The 
maximum allowed value of <f> therefore corresponds to 
E p = S (2A/J) 1/2 SA - 2SK 1 _ sin 2 = 0. This is related 
to the current through the terms in K± £ i Si y 2 which 
are linear in b i and b\. These are imaginary and can be 



3 



written as -i(J /2A) 1 / 4 K ± (sin2fa(bl-b w ). This cancels 
H t j if (hja 2 /2eS) = (5/2) (J/2 A) 1/2 K± sin 20 whence 
eliminating </>, assuming this is small, gives a critical cur- 
rent 



and gives 



Jcl 



2eS 2 
ha? 



J 

2A 



1/4 



K ± SA 



(6) 



If K± is small there is an intrinsic maximum current, 
j r = {eS 2 /ha 2 ) (J/2A) 1/2 K ± , for <j> = tt/4, which will 
be reached before E p = 0. It follows for large SA that 
jr < jci and K± is negligible. In order to conserve an- 
gular momentum, the wall must now produce magnons, 
i.e., be time dependent. Sought again are current carry- 
ing eigenstates. To the j = wall is added a twist char- 
acterized by k(z) = d(j)/dz and the solution is rendered 
stationary by a rotation of angular frequency ui about the 
new axes 9i and fa. The final axes of quantization then 
correspond to Oi = 9i+fa and fa. The effect of such a ro- 
tation is to add an effective magnetic field hw/ g^ig which 
has an component sin flihuj/ gnu perpendicular to the 8j, 
fa axis. The twist of the bond i — > i + 1 by (d<fii/dz)a 
is introduced by i? M+1 = eWW'Si'^i s *'*. Equilib- 
rium, in the rotating frame, implies zero coefficients of 
both b\ and b i . The easiest fashion by which to establish 
these equilibrium conditions is to observe the local Sf x = 
(2Sf/ 2 (l/2)(b t + b\) while Sf y = (2S) 1 /\l/2t)(b i - b\). 

Given that, e.g., R y = e tSi y se rotates the ith spin by a 
small 69 it is easy to show by evaluating [Sf y ,'H] that the 
coefficient of Sf x in H must be (l/2S)(dE / d9 t ) where 
E{9i, fa} is the already defined classical energy. The 
exchange part of this is Ej = — J^Jcos 6i cos Gj+i + 
sin6i sin6i + i(cos0i cos0 i+ i + sin^j sin^+i)]. A small 
rotation by 5<j> is generated rather by R z = e tSi ' Scl> 
where S lz = cos 0^ + sin 6^. However [S{ X ,H] 
contains no c-numbers and so this rotation determines 
the coefficient of Sf , i.e., this is (1/2S sin 9i)(dE/d fa). 
That the coefficient of S[ x be zero implies [AS sin20i + 
JSa 2 (d 2 @i/dz 2 ) + (JSa 2 )(dfa/dz) 2 sin29 4 + fiwsin/Sj] 
= 0. This can be re-integrated by 0j to give back the 
energy to within a constant, i.e., it is implied that: 



[(AS- 



JSa' 



fc(z) 2 )sin 2 9 J + (1/2) JSa 2 (d@ i /dz)' 2 



-hw(cosf3i - 1)] = 0. (7) 
Using dEj / d<f> for the coefficient of Sf is obtained 



dz 



JSa 2 



(sin Q^^-kiz) sin 2 6 4 
oz 



0. 



Integrating this with respect to z results in an equation 
relating currents: 

(J5a)fc(z)sin 2 e 4 + ^cose j = -^ (8) 
eS eS 



k(z) = 



hja 2 



2(JSa)cos 2 ^ eS 



(9) 



which, with Eqn. (7), gives an equation for Oi. How- 
ever, except to the extreme right of the wall, Qi s=s 6i 
and Eqn. (7) gives directly while Eqn. (9) determines 
the twist k(z). For the extreme right Eqn. (9) reduces 
to (JSa)k s f3 2 = (hja 2 /eS) which equates spin currents 
and ultimately determines amplitude of the magnons 
leaving the system to the right, while Eqn. (7) reduces 
to hw — 2 AS + JSa 2 k 2 which, correctly, equates the 
frequency of the rotating frame to the magnon energy. 
The whole solution is then parameterized by the value 
of k s . This, in turn, is determined by requiring that the 
absolute value of the ground state energy be a minimum. 
Given that the system is large compared to the width 
of the wall, the energy to be minimized it that of the 
coherent magnons to the right of the wall. The result 
is k s a — (2A/ J) 1 / 2 which corresponds a magnon energy 
per unit site of e r = (2Aj J) l / 2 (hja 2 / e). 

A relationship between e r and the force on the wall 
is obtained by considering a wall displacement Az gen- 
erated by R(Az) = 1 + S 1 l 2 (2A/J) 1 / i (Az/a)(bl - 
b w ). For the displaced wall HR\0) = [E + 
e r (Az I a)]i?|0) , not accounting for pinning. It follows 
that AE = [H,R(Az)\ = e r (Az/a) and that there 
is force f = AE/Az = e r /a. Using [H,R(Az)} = 
S 1 ' 2 (2A/J) 1 ' i (Az/a)[H,(bl - bj] = e r (Az/a) shows 
H contains {2S)- 1 ' 2 (2A/ jy/ 4 (hja 2 /e)(bl + b w ). A 
similar term arises from the linear part P = 
SSAsin9o cos9o(2S) 1 ^ 2 (bl ) + b ) of the pinning term, i.e., 
[P,6J = S 3 / 2 (J/2A) 1 / 4 S6Asm9 cos9 (d9 /dz)a. Cor- 
responding to 9 = 55°, the maximum pinning force, is 
(4/3\/3)S' 3/2 (5A(^+& u ,). Finally, the magnon de-pinning 
current is obtained: 



Jc2 



eS< 



3V3 V ha 2 



SA. 



(10) 



This key result is equivalent to equating the rate of cre- 
ation of linear magnon momentum to the maximum force 
produced by the pinning potential. 

The voltage across the wall is easily deduced. For large 
S, a non-polarized contact will inject majority f elec- 
trons on the left and this defines /if = m as the chemical 
potential of the left contact. For an eigenstate, the en- 
ergy to add the last 1 electron cannot depend on position 
and thus /x-f is the chemical potential for these electrons. 
One could imagine measuring /if on the right with a spin 
polarized contact. The destruction of such an electron 
involves destroying a magnon and hence /if must be ex- 
actly huj£ higher in energy than /if = /i r , i.e., the ma- 
jority and contact chemical potential to the right, i.e, 

— Mr = e V = ~ ^S for small currents. There is 
evidently a branch imbalance. 
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Current density (10 1 ^\/m 2 ) 



FIG. 1: The experimental points are taken from [3]. The 
solid line is estimated to correspond to the value of C from 
Eqn. (If) and to a j k « f.O x f0 12 A/m 2 . (This has been 
corrected for the fact that the experimental material has a 
polarization p ~ 0.7.) Evidently to within the errors the 
conduction electron angular momentum, not destroyed by the 
pinning center is, transmitted to wall motion with efficiency 
equal to p m 0.7. 

By design, Permalloy has a small intrinsic anisotropy 
and both A and K± might be expected ~ 2ttM typi- 
cally ~ 0.1T. This is several orders of magnitude smaller 
than J and it must be expected that, e.g., variations in 
the wire width lead to variations SJi and pinning. The 
largest effect occurs for the components of SJi with a 
length scale which are comparable to w the wall width. In 
this case the pinning gap E p can be as large as ~ SV AJ . 
With this it would be the case that j r < jci and the time 
dependent solution with j c i will apply. 

Having considered free motion and pinning it is easy to 
solve the problem of slow motion with a pinning poten- 
tial. If a wall is moving slowly and there are many pin- 
ning centers, the wall will adiabatically adjust to remain 
in the ground state. The wall width is w = {J /2A) x / 2 a 
and therefore the characteristic time the wall, with veloc- 
ity v, interacts with the pinning center is At — w/v and 
corresponds to an energy E c ~ h/ At ~ h{2A/ J) 1 / 2 (v / a) . 
Adiabacity requires that this be less that the energy of 
the magnons which are involved in the pinning deforma- 
tions. At the worst these have a k ~ 1 /w and the energy 
involved is ~ 2SA the gap in spin wave spectrum. Adi- 
abatically then implies v < v c ~ 2S(a/H)(AJ) 1 / 2 which 
is ~ (A/ J) 1 / 2 times the maximum spin wave velocity. 
Given the experimental values [5] correspond to v ~ 3m/s 
this criterion is surely satisfied. 

The principle of the calculation is to divide into two 
parts the rotating frames used to render the Bloch wall 
stationary, one of which corresponds to the translational 
motion and a second which reflects the pinning. The dis- 
placement corresponds to the rotation operator R(Az) = 
rii(l + {l/2)(d6 l /dz)Az{2S) 1 / 2 {b\ - bA) where Az = vt. 
Given the adiabatic approximation is well satisfied the 
rotations R(vt) can be absorbed into slowly time depen- 
dent conduction electron operators R(vt)c i R^ 1 (vt). In 



this moving frame, the effective current density is reduced 
to f = j(l — -^). Then following the calculation for 
pinning, the remaining rotations R? are those required 
for the pinned ground state. The resulting j' is given 
by Eqn. (10) and is determined by the pinning strength. 
The difference between the actual current and that which 
can be pinned is j — j' , and this in turn, determines v 
and its relationship with the current j. The velocity v is 
given by: 

v = C(j-j k ); C=^ (11) 

where, as stated, jk is the j c given by Eqn. (10) except 
that average rather than maximum pinning strength is 
involved. The important result here is that C is inde- 
pendent of jk, i.e., that part of the angular momentum 
current which is in excess of that converted to magnons 
by the pinning potential is 100% converted into motion 
of the wall. In Fig. (1) this prediction, corrected for 
p = 0.7, is compared with the experiments of Yamaguchi 
et al. Experimentally it is observed that there is a min- 
imum velocity of ~ 3m/s for which uniform motion is 
observed. It is reported that no motion of the wall is 
seen for currents smaller that 1.0 x 10 12 A/m 2 and taking 
this to be jk implies the solid line shown in the figure 
with a gradient which would correspond to ~ pC. What 
is strange is that the same velocity is reported for a cur- 
rents in the range 1.15 x 10 12 A/m 2 to 1.25 x 10 12 A/m 2 . 
It is unfortunate that it is not possible to follow the sys- 
tem to higher velocities in order to see if the initial points 
have to do with threshold effects, e.g., reflecting the first 
order nature of the transition to the magnon producing 
pinned state. The wall will be displaced during the time 
that the magnon amplitude is developed and, in a pulse 
experiment, this might be interpreted as depinning. 
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